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ABSTRACT 

One of the most important and uncertain stages in the binary evolution is the common envelope (CE) event. 
Significant attention has been devoted in the literature so far to the energy balance during the CE event, expected 
to determine the outcome. However this question is intrinsically coupled with the problem of what is left from the 
donor star after the CE and its immediate evolution. In this paper we argue that an important stage has been 
overlooked: post-CE remnant thermal readjustment phase. We propose a methodology for unambiguously defining 
the post-CE remnant mass after it has been thermally readjusted, namely by calling the core boundary the radius in the 
hydrogen shell corresponding to the local maximum of the sonic velocity. We argue that the important consequences 
of the thermal readjustment phase are: (i) a change in the energy budget requirement for the CE binaries and (ii) a 
companion spin- up and chemical enrichment, as a result of the mass transfer that occurs during the remnant thermal 
readjustment (TR). More CE binaries are expected to merge. If the companion is a neutron star, it will be mildly 
recycled during the TR phase. The mass transfer during the TR phase is much stronger than the accretion rate during 
the common envelope, and therefore satisfies the condition for a hypercritical accretion better. We also argue that the 
TR phase is responsible for a production of mildly recycled pulsars in double neutron stars. 

Subject headings: binaries: close — stars: evolution — X-rays: binaries — pulsars: general 



1. UNCERTAINTY IN THE COMMON ENVELOPE THEORY 

In the standard treatment of commo n envelope (CE ) 
outcomes via the "energy formalism" () Webbinld Il98l ■ 
the final separation of the binary is determined by equat- 
ing the binding energy of the (shunned) envelope i^bind 
to the decrease in the orbital energy Eo^h'- 
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Here ai and at are the initial and final binary separations, 
mi and TO2 are the initial star masses and nic is the final 
mass of the star that lost its envelope. 

-£'bind is considered to be the sum of the potential en- 
ergy of the envelope and its internal energy, and can 
be found directly from stellar structure for any accepted 
core mass (there are also modifications for -Ebind, where 
ionization energy or e nhanced w inds are tak en into ac- 
count, e.g. lHan et alH TQQS. 2002; ISokeH [2004): 
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Here A is a parameter introduced to fit -Ebind! it charac- 
terises the donor envelope central concentration, rric is 
the mass of the removed giant envelope and is commonly 



assumed to be rric 



mi 



Ri is the radius of the 



giant star at the onset of CE, and e is the sum of the 
specific internal and potential energies. 

For the final balance of energy, one more parameter 
is introduced, acE, to measure the energy transfer effi- 
ciency from the orbital energy into envelope expansion: 



/ Gmcm2 Gmim2\ Graim. 



V 2af 



2ai 



R^ 



(3) 



1 University of Alberta, Dept. of Physics, 11322-; 
monton, AB, T6G 2E7, Canada 



Ave, Ed- 



We anticipate that introduction of the two parameters 
introduced accordingly two uncertainties. It is common 
to remove these uncertainties at the same time, consid- 
ering the product of acE and A, by means of compari- 
son of observations with the binary population synthe- 
sis calculations, where the product of the two parame- 
ters is varied to match the observations. However, this 
approach has shown inconsistencies with the observa- 
tions, especial ly large for the formation rates of black 
hole LMXBs f|Podsiadlowski et al.l[200l lJustham et al.l 
120061). In particular, fo r LMXBs this required acEA > 2 
(lYunge lson et al."2006'). although in massive giants A <^ 
0.1 (jPodsiadlowski et al. 2003), and acE is bound to be 
< 1. 

The other way to reduce uncertainties is to consider 
them separately, e.g. one can try to determine an 'ac- 
curate' value of A from stellar structure calculations. 
It is then crucial to be precise about the definition of 
the core - should only the hydrogen envelope be re- 
mo ved, or together with the H-burning shell, and so 
on (|Tauris fc Dewil 120011) . Without knowing what ex- 
actly counts as the core and which material ought to be 
ejected, the inferred A can vary by a factor of several 
from this uncertainty alone. 

The physical reason for this variation is that in gi- 
ants, within the hydrogen shell, the potential is strongly 
increasing towards the core. The uncertainty increases 
as the mass of the donor increases, and changes from a 
about a factor of 2 in intermediate mass stars at early 
giant stage to a factor of 20 and m ore for well-evolved 
massive stars (|Tauris fc Dewil [200l . 

We stress that neither observations nor theory provide 
now a strong constraint on what post-CE remnant mass 
should be at the moment when the dynamical phase of 
the ejection ends. It does not have to be the same as the 
mass of the remnant that we observe now, e.g., in double 
white-dwarf (WD) systems or in sub-dwarf B stars: some 
remaining post-CE hydrogen-rich material can be easily 
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removed through strong winds similar to those on hor- 
izontal branch, or asymptotic giant branch, or in Wolf- 
Rayet stars etc. Between the dynamical phase and long- 
term evolution, the core will readjust itself on a thermal 
time-scale, and this has not been addressed. Here, we ad- 
dress the problem of what the post-ejection mass could 
be, different regimes in which a post-CE remnant can 
shed its remaining hydrogen-rich mass and the conse- 
quences for a companion due to post-CE mass transfer. 

2. THE POST-EJECTION REMNANT 

2.1. The divergence point 

In SPH simulations of physical collisions between a RG 
and a neutron star (NS) it has been found that not all 
hydrogen material is ejected along with the envelope - a 
tiny layer of hydrogen, from the H-burning shell, remains 
([Lombardi et al.l 120061) . This event is not directly com- 
parable to a typical CE event in a binary as at the time 
of initial approach, at periastron, the H-burning shell of 
the donor could have been in the immediate Roche lobe 
of the intruder. This magnifies the mass-loss from the 
H-burning shell and as such can decrease the mass of the 
post-CE remnant compared to a typical CE, where this 
shell might never be in the Roche lobe of the spiraling- 
in companion. This example makes clear that even in 
a dynamical CE some hydrogen-rich material always re- 
mains. 

On the other hand, studies of the evolution of stripped 
cores of low-mass RGs, have shown that there is a min- 
imum "envelope" mass Smc^min that has to be left on 
the core in order for the star to reexpand; if less mass 
is left on the core the star will contract and become a 
WD (jD einzer &: von SengbuschI [19701 ). This expansion 
or contraction of the remaining shell occurs on the ther- 
mal timescale of remaining layer, rth- 

It is plausible therefore to suppose that there is a 
unique "divergence" point nid inside the hydrogen burn- 
ing shell, such that if a post-CE star has any mass above 
this point, the star will continue to expand on rth- If its 
final mass is less than md, the star will shrink, also on 
its Tth- We recognise that rth might mean different val- 
ues in the case of degenerate core (applicable only to the 
remaining shell) or non-degenerate core (where it likely 
to depend on the core conditions). We expect that the 
material above the divergence point, if left, will expand, 
in order to obtain thermal equilibrium, but is not re- 
quired to escape to infinity without an additional energy 
source (such as the orbital energy). During this thermal 
readjustment (TR), it may also fill its Roche lobe. 

2.2. Calculations 

We tested this idea of "divergence" point on giants 
of several initial masses (1, 2, 10, 20, 30 Mq). The stars 
were evolve d using the stella. r code and input physics 
described in llvanova fc TaamI ()2004f ). This code is ca- 
pable of performing both hydrostatic and hydrodynamic 
stellar evolution calculations. For a Roche lobe over- 
flow evolution in binaries, it finds mass loss rates im- 
plicitly. Massive stars, where wind loss are important, 
were e volved with wind loss rates according to lVink et al.l 
(|200lD, or, where Vink rate s are n ot applicable, accord- 
ing to iKudritzki fc ReimersI (|1978D . 

For each initial mass, we chose 2-4 evolutionary states 
within the giant stage with different hydrogen-exhausted 
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Figure 1. Comparison of the local thermal time-scale Ttii(n^) and 
the local dynamical time-scale T^yn{m) with the mass-loss time- 
scale TTviL- Shown in the case of 18.5A/0 (ZAMS mass 20Mq, 
considered when R = 750 Rq). 

core masses mx . As during the advanced evolution stages 
stars can shrink, we ensured that the chosen giants had 
expanded to their current radius for the first time. On 
these giants, we imposed very fast ( "adiabatic" ) mass 
loss, IMQ/year. Such timescale for the mass-loss tml - 
about several initial binary orbits - is comparable with 
a fast CE event. The lower bound on tml ?J lyr should 
be clear as a CE event has to happen over at least one 
binary period at the initial Roche lobe overfiow. 

We do not imply that a CE ejection features a con- 
stant fast mass loss, and also do not study the reaction 
of the outer (convective) envelope. We are interested in 
the reaction of the inner layers, which are most likely to 
remain after the envelope ejection has occurred. 

For each mass coordinate, let us compare tml with the 
local thermal timescale rTH(™) = £^bind(w)/L(?7i) and 
the local dynamical timescale rdyn('7i): The mass-loss 
and the star evolution will be adiabatic if tml ^ ''th("i). 
The evolution can be described by hydrostatic approxi- 
mation if Tml ^ Tdynirn), as a star will always acquire 
its hydrostatic equilibrium within a dynamical time, and 
its state at the hydrostatic equilibrium is defined by its 
thermal structure. 

For most stars, tml is much shorter than any local 
thermal timescale (see Fig. [1}. We note however that 
in the inner layers that are close to the cores of our 
most massive stars (20 and 30 Mq), the complete mass- 
loss sequence can take up to 10% of the local thermal 
timescale of a few hundred years. Thus, even such a fast 
mass-loss produces only approximately adiabatic evolu- 
tion: some thermal evolution proceeds and is expected to 
be responsible, in particular, for some expansion of inner 
non-degenerate layers during the CE phase. As a san- 
ity check, we calculated additional mass-loss sequences 
for massive stars, with faster and slower mass-loss rates, 
and found only minor differences in the region of interest 
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between the runs with 0.1, 1 and 10 Mqyt'^. 

On the other hand, local dynamical time-scales are 
longest at the surface and significantly shorter for in- 
nermost layers (Fig. [1]). Tdyn(TO) is comparable by the 
order of magnitude to tml in the outer layer of the mas- 
sive giants. Tdyn{fn) is however by 3 or more orders of 
magnitude smaller than tml in the Helium rich layers, 
and closer to Hydrogen exhausted core, it is ~ 10~^rML, 
even in our most massive considered stars. Therefore, 
although the evolution of the outer layers is indeed de- 
pendent on the inclusion of hydrodynamical terms, the 
inner layers always have enough time to regain hydro- 
static equilibrium, and are therefore insensitive under 
the adopted mass-loss rate. In summary, we find that 
for studies of the thermal reaction of inner layers that 
will form the remnant after the fast envelope ejection, 
the hydrostatic version of the code is sufficient. 

As a result of mass loss evolution, we obtained se- 
quences of (post-CE) remnants with different final (post- 
CE) masses, each of which then was evolved for several 
Tth, to check if this post-CE star is expanding or con- 
tracting. We note that in our code the value of post-CE 
mass could be resolved no better than pre-CE resolu- 
tion in hydrogen shell, this is specifically important for 
low- mass giants (e.g., we have about 50 mesh points per 
O.O2M0 H-sheU in 2 Mq RG with a mx = 0.52Mo). 

To summarize, we separate the CE event into two 
stages: one resulting in the envelope ejection, and the 
subsequent thermal readjustment of the remnant. The 
latter part is a distinct phase unless the spiral-in (includ- 
ing the ejection of the envelope) takes place on a time- 
scale comparable to the shortest thermal time scale, sev- 
eral hundred years; and has not been heretofore treated 
in the literature. 

Finally, an admonishment is in order: several esti- 
mates exist for the CE duration, neither one of them can 
boast conclusive observational evidence or indeed self- 
consistency. E.g., a 'slow' CE c ould last for 100 years 
and longer (jPodsiadlowskil 120011) . We can not justify 
which CE evolution timescale is more appropriate, and 
this is not the purpose of this paper. We concentrate 
on the 'fast' event, however, we see no reason why our 
results should not be applicable in the 'slow' case: the 
core reaction will be similar, albeit lagging by the time 
the ejection takes. We also note that a 10-times slower 
loss rate did not produce a significant difference in our 
calculations. 

2.3. Degenerate cores 

Indeed, as in previous studies, we found that every 
low-mass giant with a degenerate core has a unique di- 
vergence point mdiv such that if post-CE mass is less 
than TOdiv, it contracts on Tth. AH post-CE remnant with 
masses above mdiv expand, create new outer convective 
zone and keep expanding even after Tth. 

After locating mdiv, we analyzed pre-CE giants struc- 
ture to find what characteristics these points had in ini- 
tial giants, before the stripping began. For this, all gi- 
ants, including massive, were used. We noticed that 
among all the giants, mdiv could have initially a wide 
range of hydrogen content, X = 0.08 — 0.58, and so a 
criterion involving specific constant value of hydrogen 
abundance could not be satisfactory. Similarly, another 
criterion discussed in the literature - the location where 



the e nergy generation rate is maximum ()Tauris fc Dewil 
I2001D — does not coincide with the divergence point. 
We found that in the considered models mdiv is close 
to the "maximal compression point" mcp, which is the 
mass zone with the maximum value of 'compression' P/ p 
in the hydrogen shell. 

2.4. Non-degenerate cores 

For massive giants, as previously, a post-CE remnant 
also has divergence point that corresponds to the min- 
imum post-CE expansion of the remnant, although the 
overall response is different from the case of the giants 
with degenerate cores: 

• For all possible remnant masses with m < mcp , the 
core slightly adiabatically expands during the fast 
adiabatic mass loss. Once we stop the mass loss, it 
can very slightly (a few per cent) expand and then 
shrink dramatically, becoming smaller than it was 
before the CE. 

• For larger remnant masses, as previously, the con- 
vective envelope is re-formed, and the star remains 
as an extended giant for a while. The envelope can 
become larger than the pre-CE gianfl 

• For the intermediate range of remnant masses, 
above the mcp, but below the boundary where the 
convective zone redevelops, the core experiences a 
pulse on ~ Tth, being able to expand by up to few 
hundred times more than this mass had as a ra- 
dius coordinate before the CE. After the pulse, the 
post-CE star shrinks significantly, also becoming 
smaller than prior the CE. 

To illustrate these three types of the response, in Fig. [2] 
we show the typical case of 9.75Mq (ZAMS mass IOMq) 
star, taken when it had radius of 300 Rq. For compar- 
ison, we show 18.5M0 (ZAMS mass 20Mo) star with 
radius of 750 Rq (see Fig. |3]). Even though this giant 
has profile of hydrogen qualitatively different from the 
considered above IOMq star, it shows similar behavior. 
The main difference with a 10 Mq star is that there is a 
more contrasting response between inside mcp and out- 
side it. In a 3OM0 giant this difference even stronger as 
mcp is located just below the bottom of hydrogen burn- 
ing convective zone. 

2.5. The adiabatic response 

We recognize that the response we discuss above is 
non-adiabatic, but is the equilibrium response (the one 
that a star experiences in order to obtain its thermal 
equilibrium). Another important response to consider is 
the reaction of a star when the dynamical event ends, 
true adiabatic response. It is established that adiabatic 
response of the surface lay ers depend on whether they 
are radiative or conve ctive (jHiellming fc WebbinkI [TQSTI : 
iSoberman et al.|[l997[ ). In particular, radiative layers on 
dynamical timescale tend to shrink, and convective layer 

Here, we can not fully separate the post-CE TR expansion 
from a normal stellar evolution along a giant branch, however we 
find that a post-CE star obtains after the CE a large radius faster 
than it would have otherwise. 
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Figure 2. Size of a post-CE remnant for a 9.75Mq star (ZAMS 
mass 10A/q, considered wlien R = 300 Rq). i?pro— CE is tlie ra- 
dius coordinate of each considered remnant mass between tlie mass 
loss, JJpost— CE is the radius of a post-CE star when the fast adi- 
abatic mass loss stopped and i?max is the maximum radius that 
this post-CE star had obtained within (time to reach maxi- 
mum ranged from 100 to 2,100 years). Shown are the ratios of 
CE and Rpre-CE (blue line, open circles correspond to each 
calculated model), i?max and Rpre— CE (green line, solid circles), 
Rmux and Rpost— CE (red line, open squares). X = is the mass of 
the hydrogen exhausted core, £n = max is where the nuclear burn- 
ing has maximum energy generation rate, rrtcp is the "compression 
point" (where P/p has the maximum value within the layer be- 
tween X = and rribcz) and rribcz is the mass coordinate of the 
bottom of the convective zone in the pre-CE star. 
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Figure 3. Size of a post-CE remnant for a 18.5Mq (ZAMS mass 
20 , considered when R = 750 Rq ) . Notations as in Figure [2] 

remain the same or expand. We checked for all consid- 
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Figure 4. The initial binding energy (solid red line) and available 
orbital energies in the initial and final configurations (dotted green 
and dashed blue lines) for an 18.5 A/q star (ZAMS mass 20 Mq) 
with R = 750i?Q. When calculating the post-CE AEaih^ the sizes 
of the post-CE remnants for every given mass were used to define 
Of. The mass coordinate is normalized to rrtcp (in this star rricp = 
6.47M©). The AE^rh are calculated assuming that the companion 
mass is 1.82 Mq. 

ered model the location of the divergence point and found 
that all of them are located within initially (pre-CE) ra- 
diative layers. We conclude that immediate response for 
mass removal to the divergence point is always a shrink- 
age and therefore does not affect our conclusions based 
on the thermal response. 

3. CONSEQUENCES FOR THE ENERGY BUDGET 

As the core expands during semi-adiabatic mass-loss, 
a surviving binary must be wider when core expansion is 
taken into account. Thus more b inaries will merge (thi s 
result is opposite to the claim in'De love fc Taamll2010( ). 
As an example, the giant in Fig. [3] will easily survive a 
CE event with a NS of 1.4M0 (with acE = 1) if its core 
did not expand. Setting core mass < TAAMq will satisfy 
the energy budget to create a compact binary. However, 
if one takes into account the core expansion, the binary 
will merge: for all core masses above TOcp, the remnant 
will expand significantly and overfill its Roche lobe. A 
minimum companion mass, for which both energy budget 
and post-CE core size are taken into account, will be 
^ 1.82M0 and the giant core in this case should have 
been removed to at least rucp (see Fig. |4]). It can be 
seen from this Figure that for a fixed companion mass, 
there is a unique solution where available orbital energy 
can exceed binding energy if the remnant expanded after 
mass loss, whereas a non-expanded remnant gives a wide 
range of the possible core masses, from 5.95 to 7.5 Mq. 

Let us introduce "the energy expense" , the difference 
between the required energy £'bind("T^) and the available 
orbital energy (5i?orb, normalized per £'bind('™)- 

r _ £^bind("^) ~ (£^orb,i " £^orb,f) , .s 



Mass and fate of the post-CE remnant 



5 



1 - 




< 

\ 

-1 - ^3=1 Mg 

II 



-2 - 



0.8 1 1.2 1.4 

Figure 5. The energy expense in the the layers that can become 
a post-CE remnant for a 4.75 Mq star (ZAMS mass 10 Mq) with 
R = 400i?Q. The two curves are for different companion masses 
(as marked). The mass coordinate is normaUzed to mcp. 
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Figure 6. The energy expense in several giants with initial mass 
20AfQ. The mass coordinate is normalized to rricp. Energy ex- 
penses are calculated using different assumed companion masses, 
to bring them to the same value for the energy minima. 



It is the (normahzed) excess energy available to the enve- 
lope after all the matter above the given mass coordinate 
has been removed. Of course, a positive signals that 
the removal process is not possible from the energy con- 
siderations. Fig. [5] shows distribution of as a function 
of the mass coordinate. For these calculations, E'orb.f 
was assumed to be at the Roche lobe limited orbit for 
the post-CE remnant (semi-adiabatic expansion is taken 
into account), so that available orbital energy is at its 
maximum. 

Note that mcp the energetically optimal position to re- 
move the envelope to: it is the equilibrium point of the 
generalized force 9m(^'bind+^'orb)- The shape and the lo- 
cation of the minimum of only depends on the donor's 
energy profile up to the Roche radius, and does not de- 
pend on the companion mass. The latter only determines 
the magnitude of the energy excess (see Fig. [5]). 

Thus it is easy to determine the minimum mass of a 
companion which allows the survival of the binary: as it 
can be seen from the Fig. [5l it is such that will result in 
removal of mass to Wcp, precisely. 

We verified that the coincidence of the minimum of 
the energy expense with shedding the envelope to about 
TOcp for minimum likely companion mass holds for many 
of the studied giants, though does not hold for giants 
that are early on the giant branch which only recently 
develop convective envelopes. E.g. in an early 10 Mq we 
observe one more energy minimum, at a higher core mass 
~ 2.45M0; the local energy minimum at rricp nonetheless 
holds. In more massive early giants, the energy expense 
minimum is between rricp and nix (Fig.[6]), its location is 
closer to the location of e„ = max than to TOcp! the star 
still has the same reaction on expansion or contraction 
with respect to other stars. 

4. POST-CE MASS TRANSFER 



4.1. Consequences for the final post-CE mass 

Let us consider what happens if the companion was 
massive enough so that during CE not all mass to mcp 
had to be removed. In this case a CE would end with 
a binary separation such that Roche Lobe overflow for a 
post-CE remnant during its TR will follow. To consider 
this, we took a post-CE model showed in Fig. |31 con- 
sidering its remnant with mpost-CE = 6.84AfQ (larger 
than its rricp = 6.47 Mq). This post-CE remnant is then 
placed in a contact binary with an arbitrary companion 
of 3 Mq: this mass self-consistently satisfies the energy 
required to shed the mass above 6.84Mq. For the mass 
transfer (MT) we can choose a fully conservative or fully 
non-conservative mode 0. 

As expected, the post-CE remnant rapidly expanded 
and started the MT; initially at a very high rate (~ 

1 — 5 X lO^^AfQ/yr), in accordance to rth of the re- 
maining hydrogen rich layer. After removing most of 
the layer above rricp, it slowed down to rth of the core 
(- lO^^Mo/yr). The MT continued till MT rates be- 
come comparable to the TR time-scale so that the core 
can shrink faster than it expands due to mass loss. In 
this particular example the core reached almost exactly 
the divergence point, shedding ~ O.4M0 during the MT 
so that the final mass was rricp- We also performed a MT 
calculation in the fully non-conservative regime. In this 
case the final mass of the core after the TR phase is the 
same as in the conservative calculations. 

We also considered the case when the companion is 
a NS. We note that with a 2OM0 donor then, even for 
acE = 1, the energy requirements for envelope ejection 
would not be satisfied if the core expands as much as 

^ Partial conservation, limited to Eddington rates can be con- 
sidered as well, but with the MT rates that we find and describe, 
this case does not differ much with the fully non-conservative case. 
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we find after our fast mass loss (i.e. the binary would 
merge). The final mass of the remnant of the massive 
giant is the same, nicp, as for the 3 Mq companion. If the 
mass transfer is fully conservative, the NS is presumably 
spun-up, as it would accumulate 0.34AfQ. Again, in the 
case of a fully non- conservative regime, we find that the 
final mass of the post-CE remnant is mcp. 

Next we considered a system with a 10 Mq giant (same 
as shown in Fig. [2]), considering as the post-CE core 
2.54 Mq (mcp = 2.04). In our MT simulations with a 
NS companion, less material above mcp has been trans- 
ferred, only 0.24 Mq, the final remnant mass is 2.3 Mq. 
It might be connected to the fact that 10 Mq star does 
not have such a sharp profile as a 20Mo in the post-CE 
thermal pulse zone, and its post-CE expansion is more 
fiatter until about this mass (see Fig.[2|). With a smaller 
companion mass, more of the post-CE remnant mass is 
stripped off. 

4.2. Consequences for the companion 

The MT rates that we encounter in the post-CE TR 
phase are highly super-Eddington and we face the obvi- 
ous question whether the MT is approximately conser- 
vative or almost non-conservative, as this is crucial for 
a companion. The question what happens if the mass- 
accretion rate on a NS exceeds Eddington limit has been 
discussed extensively in the literature, in particular, the 
regime in which it exceeds A^Edd by many orders of mag- 
n itude. 

iBegelmanI (|1979f ) showed that if the accretion rate is ex- 
tremely high, few times 1Q~'^Mq yr~^ , then within some 
volum e (the "trapping radius", e.g. iKing fc Begelmaiil 
|1999() around the star the diffusion of photons outward 
cannot overcome the advection of photons inward. While 
a black hole can swallow all the material in this case, if 
the accretor is a NS, radiation pressure near the NS's sur- 
face resists inflow in excess of the Eddington limit, likely 
leadin g to creation of a Thorne-Zytkov object. Blondini 
(|1986D has also found that when MT rates exceed the 
Eddington rate by 10^ x L-^dd/c^ or more, the accretion 
proceeds in a hypercritical regime. 

Hypercritical accretion was then argued to be respon- 
sible for such efficient material accumulation during a 
CE event, that a NS is likely to convert to a black hole 
(IChevalierill9"8l . IBrownl (11991 used this argument to 
understand double NS formation. He showed that in- 
deed in a CE event the Bondi-Hoyle-Lyttleton accretion 
rate is about 10'' x MEdd and a NS can accumulate up 
to 1 Mq . He argued that in this case, considering that a 
number of the discovered double NS have masses closer 
to the lowest possible NS mass limit, a double NS can 
be formed only from a binary with almost similar initial 
masses, evolving then via double CE event, before either 
o f the NSs was forme d. 

IHouck fc Chevalieil (jl991[ ) considered neutrino losses 
during accretion on a NS. They studied in detail the 
regimes of the mass accretion 10"'' ^ < IC* Moyr"'^, 
and found that radiation diffusion becomes important 
when the accretion rate falls below lO~^M0yr~^, for 
smaller rates the radiation pressure can not support an 
envelope around NS surface and can not cease the in- 
fall of the material. We note that the accretion rate 
that separate the hypercritical accretion with the accre- 



tion when the radiation diffusion dominates in this case 
(lO^^Af0yr^^) is higher than the one found to work dur- 
ing a CE event (2 x lO-^MQyr-i). 

If the latter estimate is more proper than in the stud- 
ies listed above, then it is possible that a CE hyper- 
accretion, as having too low mass accretion rate, does 
not lead to a significant accumulation of the material. 
Post-CE core expansion, however, in either case can lead 
to a hyper-accretion regime, as during this thermal pulse 
it provides a much higher mass accretion rate. This can 
lead to a NS spun up. In our calculations, MT rates 
exceeded lO~'^iW'0yr~^ long enough to accrete in hyper- 
critical regime on a NS O.29M0 in the case of 2OM0 giant 
and O.O9M0 in the case of a lOM©. 

We note that the observed double NSs are ge nerally 
mildl y recycled, having periods 0.024—2.7 seconds jStairsI 
I200I . The mass distribution of those with mass mea- 
surement errors < O.O2M0 are such that the difference 
bet ween the masses of the NSs is < O .1M0 (e.g., see data 
in lStairsll2004t iKiziltan et al.ll2010D . Their location on 
the P — P diagram for galactic field NSs is also interme- 
diate between the millisecond p ulsars and non-recycled 
ones (jArzoumanian et al]|1999f) . and the post-accretion 
perio d is likely to be not millisecond (jLorimer et al.l 
^200^. It could be a sign that the very rapid mass trans- 
fer followed the CE and preceding the second NS forma- 
tion is not capable to fully spin up and efficiently reduce 
the magnetic field on a NS that was formed first. Same 
mechanism can lead to a formation of mildly recycled bi- 
nary pulsars with low-mass white dwarf companions (|Lil 
,2002; DclovG .20081. 

5. CONCLUSIONS 

We analyzed a set of giant models with respect to their 
likely post-CE response. Although our set was not ex- 
haustive, it did exhibit a clear trend that allowed us to 
conclude that: i) every giant has a well-defined post- 
CE remnant after it has been thermally readjusted, most 
likely given by the divergence point (see discussion be- 
low); ii) the divergence points, at the current resolution, 
are best approximated by the point in the hydrogen burn- 
ing shell that had maximal compression (local sonic ve- 
locity) mcp prior to CE. This definition allows us to find 
quickly a post-CE core mass for any giant without per- 
forming mass loss calculations. 

We remark that this divergence point does not nec- 
essarily mark the final mass of the remnant (e.g., the 
stellar wind in He rich stars could quickly and effectively 
remove the remaining hydrogen-rich envelope) or imme- 
diate post-ejection mass, however it marks the mass after 
the thermal core readjustment. 

The post-TR core, defined by mcp, most likely co- 
incides with the post-ejection mass in low-mass giants, 
where reestablishing of a convective envelope for masses 
above the divergence point happens on a few dynami- 
cal timescales, leading to another (likely unstable) MT 
event. As the -Ebind of remaining shells during this pe- 
riod has not been changed much compared to the pre-CE 
state, the whole sequence of events can be considered as 
one CE event with a final core being mcp. 

For giants with non-degenerate cores the situation is 
more complicated. It is not possible to claim that the 
divergence point defines the post-CE core (dynamical 
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phase) uniquely from the energetic budget point of view. 
However, the divergence point does appear to define the 
remaining post-TR core, if the post-CE configuration al- 
lows Roche lobe overflow for the post-CE remnant during 
the core TR. During this period, (stable) MT can pro- 
ceed, resulting in the enrichment of a companion star 
with material from the hydrogen-burning shell. 

In all cases, the post-CE remnant of a giant with a 
non-degenerate core has expanded by few times by the 
end of TR. Remnant masses greater than lead to 
greater expansion; for each particular giant star, the min- 
imum possible size change is for a remnant mass rricp. 
The expansion of the remnant means that less orbital 
energy is available to eject the envelope than if there was 
no expansion, since the surviving binary must be corre- 
spondingly wider. This leads to a reduction in the num- 
ber of binaries which can survive CE. So fast CE allows 
more binaries to survive the end of the dynamical phases. 
Slower, self-regulating CE leads to more mergers. 

We note that for both types of giants rricp firmly repre- 
sents only a maximum post-TR core mass, as we can not 
fully rule out that the dynamical phase will not already 
have removed mass below mcp. However, we argue that 
such extra mass loss is not likely to happen if, during 
the final stages of the spiral-in, the characteristic orbital 
evolution time is comparable to the core response time 
near uicp point, which is as short as 10-100 years. 

We also find that in most evolved giants, the energy 
required to shed the envelope down to Wcp is the min- 
imum energy expense: per total E'bind unit, it requires 
more orbital energy to remove either less or more of the 
mass from the expanded core. It is fully reasonable to 
remove less of the envelope (and have a bigger post-CE 
mass) once -Ebind < a-Eorb, the TR phase will then re- 
move mass down to ~ rricp. However, it is not plausible 
to remove the envelope deeper than to ~ ?7icp: for any 
remnant mass less than mcp, the difference between -Ebind 
and aEorh increases compared to their value at rricp. 

We suggest therefore that a divergence point uniquely 
defines the core in the post-TR phase, and then a slow 
wind loss phase follows the CE with almost no core evo- 
lutioifl A companion can be spun up during the MT ef- 
fected by the remnant's TR. Roche lobe overfiow phase is 
short, however can lead to at least mild recycling. If the 
companion is a NS, its recycling will depend on whether 
the conservative mass transfer (due to the hypercritical 
accretion) is possible or not. We find that hypercrit- 
ical accretion is more likely during this TR phase than 
during a CE, as the mass accretion rates are significantly 
higher. It also leads to a smaller mass accumulation than 



is found to occur during a CE, corroborating the mass 
distributions of the observed double NSs. We conclude 
that the post-CE TR phase can be responsible for a for- 
mation of mildly recycled pulsars in post-CE binaries and 
specifically in double NSs. 
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* We do not include here evolution on nuclear timescale which 
will follow as usual; e.g., a core of a small enough mass can again 



become a giant, as it is customary for low-mass He stars. 



